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Abstract 

Influence of relative spin alignment on the band splitting and magnetic excita- 
tions in bilayer correlated systems is studied. Splitting occurs to be large or small 
depending on relative orientation of staggered spins of the layers. Change of the 
ground state spin configuration with doping is shown. Behavior of bilayer splitting 
in Bi2Sr2CaCu20g + s allows to suppose that superconducting transition is accom- 
panied by transformation of spin configuration of system. 

PACS: 71.10.Fd, 71.27.+a, 71.10.Hf 

One of properties of high T c superconductors - a dependence of their characteristics 
on number of the CuOi layers - shows importance of interlayer coupling in them. In re- 
cent photoemission experiments |J 0, [| the bilayer splitting of the bands and the Fermi 
surfaces (FS) in Bi 2 Sr2CaCu20s+s (BSCCO) have been observed. Unusual observation 
was a significant decrease in bilayer splitting at k ~ (ir, 0) in course of the superconduct- 
ing (SC) transition [Q. Another phenomena which is observed in the bilayer (but not 
monolayer) cuprates is the appearance in SC state of a spin resonance in inelastic neutron 
scattering j|, ||, |(| . Explanation of these phenomena is important for understanding the 
electronic structure and mechanisms of HTSC displaying many properties of doped Mott 
insulator. The main feature of latter is the Hubbard splitting of band induced by antifer- 
romagnet (AF) fluctuations. Though the long range AF order disappear at small doping 
5 C = 1 — n c ~ 0.05, the local 2D spin order probably retains in sufficiently large doping 
range including the range of superconductivity. Among arguments in favour of such local 
AF order there are data on /xSR, NMR, NQR |7], ||, large length of spin correlations 
with Q = (ir, 7r) and smooth evolution of collective magnetic excitations with doping (see 
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review jp) and others. Finally, recently a direct proof of the AF order even in SC state 
of YBCI2CU3O6.5 have been obtained || |I| in nanosecond scale. 

Detailed analysis of the AF bands and interactions of the Hubbard model have been 
done in classic work fLl] . Simple mean field method gives the overestimated values do of 
staggered spin and of the doping boundary 5 C ~ 0.45 of a range of the AF order. The band 



renormalization on base of the zero AF approximation |Tl| or the calculations of t-t'-U 



Hubbard model using the slave-boson technique [12] or the states with the valence bond 



correlations 0] - all these methods give the smaller values of do and of the boundary 
doping 5 C ~ 0.3 at which the local magnetization disappears. According [|13[ the range 
5 < 5 C includes entirely the range of superconductivity. Retaining of the Hubbard splitting 
of band at 5 < 5 C in the t-t'-U |13[ and t-t'-J |14], [15|, [16] models implies the change of the 
FS topology at the optimal doping with transition from the hole pockets to large FS. 

If such picture is true then interactions between the layers in bilayer systems may 
significantly depend on relative orientation (alignment) of local staggered spins of the 
layers. It must be so even if the spin orientation is not stationary in strong sense, but 
is retained only on some time interval. In case of small difference in the energies of 
various configuration of spin alignment one can expect a strong influence of doping and 
temperature on the properties of magnetic excitations in system. A search of such effects 
may help to elucidate the origin of the magnetic resonance in neutron scattering in YBCO, 
BSCCO. 

Present work summarizes the results of variational calculations of bands, full ener- 
gies, phase curves T c (5) and the band splittings in the bilayer t-t'-U Hubbard model. 
Calculations are made on base of the correlated states with the valence bond (VB) for- 
mation - the band analog of the RVB state of Anderson. Previously ]TB] it was shown 
that the correlated hopping interaction appearing in the effective Hamiltonian during 
the VB formation does provide the hole attraction in d- channel and corresponding d- 
superconductivity compatible with AF spin order. Here we use the same method for 
study the effects of interactions of two layers. The questions of interest are: Is the large 
doping range of local AF order retained? If so, then what are the relative alignment of 
spin systems of both layers and the bilayer splitting at various doping? 

A difference between given approach and the study ]T7|] of bilayer Hubbard model by 
the FLEX (fluctuation exchange) method must be stressed. According FLEX the hole 
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attraction is provided by long range spin fluctuations. In our approach the short range 
AF correlations are decisive. In simpliest variant they correspond to the VB formation, 
i.e. formation of singlets on the bonds of the nearest centers. Unlike the averaged consid- 
eration in FLEX we make the calculations for two specific spin configurations of bilayer 
system - homogeneous states with the same or opposite alignment of staggered spins of 
both layers. We denote them as the F z or AF Z (ferro- or antiferro-magnet over z direc- 
tion) alignment correspondingly, though both cases refer to local AF spin order in each 
layer. We consider only homogeneous states, though many neutron experiments give in- 
dication on the stripe-type or spiral incommensurate structures in cuprates. Nevertheless 
the effects discovered in our calculations for two homogeneous spin structures must un- 
doubtedly be taken into account in discussing the observed bilayer splitting and magnetic 
properties of cuprates. 

We start from the Hubbard Hamiltonian of bilaeyr system H(U,t,t',t z ). The main 
intralayer interactions are determined by standard parameters U, t of Hubbard model. 
Additional interactions include the hoppings t', t z between the next nearest centers inside 
the layers and between the cites of different layers. Positive signs of standard parameters 
of strong coupling t, t', t z are defined in such way that the zero bands (at U = 0) have a 
form 

e± = —2t(cosk x + cosky) + At'cosk x cosk y =p -t z (cosk x — cosky) 2 (1) 

Expression for splitting Ae(k) = e + — e_ between the bonding and antibonding zero bands 
is derived in |l7f . 

A variational correlated state \I/ = W(a)$> with correlations of the VB type is con- 



structed |L3[ as unitary transformed uncorrelated state $ 

W(a) = exp aZ ; z = \ E ( c \na c -/ma - h.c.)(n in -„ - n im ,-a)] (2) 

u,<nm>,7 

Here 7 = 0, 1 refer to layers a,b. A choice of unitary operator H^a) with variational 
parameter a is explained in Variational treatment of problem with Hamiltonian 

H in basis of correlated state \1/ is equivalent to treatment of the effective Hamilto- 
nian H e ff(a) = W\a)HW{a) ~ H + a[H,Z] + ^-[[H,Z] in mean field approxima- 
tion. We use the most general uncorrelated states {$} of the BCS type with anoma- 
lous averages of d-symmetry and the AF spin order. The mean energy over such state 
< H >q,=< H e ff >$= H(y lu ,z) is obtained as an explicit function of a set of one- 
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particle averages {y u , z} over uncorrelated state The set of these averages contains 



the density components r 7 j =< c 7racr c 7in+ z i(J >, analogous components of staggered spin 
d 7 i = (-l) n |f| < c 7nCT c 7jn+ii(J >, anamalous averages w 7 i = (l 2 x - l 2 y )< c 7nCT c 7jn+ i ;CT > of 
d-symmetry for each layer 7 = 0, 1 and a value z = t^ 1 < T z > determining the aver- 
age interlayer hopping. For spin components we study two variants of relative alignment 
d y= i t i = Cd^7=o,h Cd = ±1- Selfconsistent procedure of the energy minimization over $ 



and of subsequent minimization over a is standard one JT3[. Details of calculations will 
be presented elsewhere. 

For main parameter we use a value U/t = 8 following from cluster derivation of single 
band model. Parameter if owing to its contribution into the band dispertion regulates the 
position Eyus — E(ir, 0) of the van-Hove singularity (VHS) in the density of state relative 
to the edge of the lower Hubbard band. In the monolayer t — if — U or t — t' — J models 
a value of t' is directly connected with the optimal doping 5 opt , since a maximum of T c 
corresponds to coincidence of the chemical potential /1 with Evhs- According the models 
[p| , p] , |15| , pfl a resonable value of 5 opt ~ 0.2-^0.25 corresponds to values t'/t ~ 0.05-j-0.1. 
These values are less than t'/t ~ 0.2 0.4 obtained from fitting of photoemission data 
and the band calculation data on base of the strong coupling scheme. Here we vary 
parameter if in limits t' = 0.05 -j- 0.3. Parameter of interlayer hopping is also varied 
as t z /t = 0.07 -h 0.3. Previous calculations without renormalization gave t z /t ~ 1/3 
[PH . Recent measurements of bilayer splitting A = 2t z = 5ek(7i~, 0) fj], @] give the value 
tz/t < 0.1. 

Fig. 1 presents typical doping dependencies of spin density do, transition temperature 
T C (S) and a difference AH = H(AF Z ) — H(F Z ) of average energies of bilayer system with 
two types of relative alignment of the layer's spins. Staggered spin density do has nonzero 
value at wide range of doping including the region of supercondutivity. Quantity AH (5) 
changes its sign at some doping. In the undoped system, 5 = 0, the configuration AF Z 
with opposite alignment of the layer's spins occurs to be lower than the F z configuration. 
At large doping S one have inverse situation up to boundary value S c , at which the local 
magnitization disappears. Negative sign of AH at 5 = is explained by a positive 
exchange constant J a b ~ 2t 2 z /U > 0. Change of sign of AH at some 5 is a result of the 
bilayer splitting of the bonding and antibonding bands, in particular, the splitting of VHS, 
and different occcupancy of them for two types of spin alignment. Maximum of AH (6) 
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corresponds to optimal doping S op t for the parent monolayer model. At this doping one 
has Ey L HS — fj, = 0. Increase of t' shifts both the 5 op t and the position of maximum of 
AH \5). Such relation is not occasional. 

According [|17] the bilayer splitting 8e£> = 2t z (cosk x — cosk y ) 2 /A of the zero bands 
has a maximum in region k ~ (71*, 0) forming the VHS in density of state. However, 
the splitting contributes to the average energy only if at this k bonding and antibonding 
bands appear to be occupied and unoccupied correspondingly. The latter takes place when 
E 1l (tt, 0) = Ey L HS = \x for the unsplit bands of the monolayer model. It remains to explain 
why a bilayer splitting manifests itself in the energy difference AH = H(AF Z ) — H(F Z ). 
The reason is in different bilayer splitting for different spin configurations. 

Fig. 2 presents the density of state of lower Hubbard band for two configurations of 
spin alignment for system with small values t', t z . Bilayer splitting is absent at AF Z 
alignment, but sharply seen in F z configuration. This is due to different behaviour of 
matrix elements of the interlayer hopping between the states of lower Hubbard band 
of each layer. At AF Z alignment these matrix elements disappear on the nesting lines. 
Different behaviour of DOS reflects itself directly on the form of the phase curves T c (5) 
for SC transition (Fig. 1). In case of AF Z or F z configuration the curve T c (5) has one or 
two maxima correspondingly. 

Thus the models with small t', t z <0.1 are characterised by 1) large T c max ~ 0.02t ~ 
116K at AF Z configuration; 2) close energies of both configurations, 3) small condensation 
energy; 4) lower F z configuration at 5 ~ 5 op t; 5) large (or zero) bilayer splitting in F z (or 
AF Z ) configuration. 

Unlike the normal state, a prediction of a lower spin configuration in the SC state is 
hardly possible. An estimate of the condensation energy and the step in a heat capacity 
occurs to be less than that observed in YBCO approximately in 5 times. For this reason 
we dicuss the spin configuration of SC state on base of the observed behaviour of bilayer 
splitting in BSCCO [1], 0]. This splitting decreases from 5e ~ 80meV in normal state 
down to 5e ~ 20meV in SC state. Such behaviour might be explained if we suppose that 
a transformation of configuration F z —> AF Z proceeds simultaneously with SC transition. 
One might suppose also that the SC transition as such is induced by the change of con- 
figuration, since such transition is accompanied by increase of the density of states at the 
Fermi level. 
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For models with large values t', t z the types of dependencies AH (5), do (5) are retained. 
Increase of t z leads to increase of \AH\ and the doping of the F z — AF Z crossover. 
Calculations confirm the wide doping range of AF local spin order for these model also. 
Value of do is greater by an order of magnitude than the value of AF spin density observed 
in SC state of YBCO || [Hj] from elastic neutron scattering. The difference, possibly, is 
connected with large distribution of the spin directions (or signs) of different biplaines in 
crystal. Significant decrease of T c and a deformation of the phase curves with increase of 
f, t z are connected with decrease of the density of state on the FS for these models. 

Influence of the spin alignment on the magnetic excitation spectrum is of interest 
in connection with discovery of the spin resonance in neutron scattering in SC state of 
cuprates and its unusual dispersion |jlQfl . A widely discussed hypothesis |L8| connects 
its origin with the so called 7r-resonance - an excitation of e — e pair with Q ~ (jr, it) 
in triplet state. But the ratio of weights of such resonance to the integral (over u) 
intensity of the non-resonance magnetic excitations would be too small. Expected order 
of magnitude of this ratio is ~ (w/do) 2 ~ 0.01. Interpretation of incommensurate patterns 
in the spin susceptibility on base of the nesting properties of FS also cannot describe a 
large intensity of spin excitations and their similarity in both the undoped dielectric 
materials and the doped metallic ones. The most probable hypothesis fl9| implies the 
common magnetic origin of spin resonance and of the incommensurate anomalies of spin 
susceptibility x"^^) i n cuprates. Both features are explained by existence of the AF 
domains or the incommensurate modulation of the local staggered magnetization. The 
stripe-type structures have been observed in LaSrCuO and they may also present in the 
bilayer cuprates. But monolayer models ||19| cannot explain the fact that the resonance 
appears only in superconducting state and only in odd channel of bilayer system. In this 
connection it is important to study the influence of spin alignment on magnetic excitations 
in bilayer systems. 

Here we study such effects using a simpliest example of homogeneous model for 
"constructing" the spin resonace. Consider the fenomenological spin Hamiltonian with 
anisotropic interlayer interaction 

H — ^2 Jo ^2 ^n^m + ^2{Jz,abS zn Sl n + J± : ab[S^ n S^, n + Sy n S yn )} (3) 
7=a,6 <nm> n 

The "^''-direction is selected by vector of local staggered spins < >= e 2 (^) 7 rfo of 
each layer (7 = 0, 1 for layers a,b) with definite alignment Q = 1 or —1, stabilized by 
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the nonspin interactions. Then the RPA consideration in the linear spin wave theory pO 



gives the following frequencies of collective spin excitations for both types of alignment 

F z : u even (q)^2Q v ^Y u u odd {q) w 2ft v /^^; (4) 
AF Z : u even {q) w 2Qy/(p + g 2 , oo odd {q) w + p i; 

Here y? = 2 + cos^ + cosg y , Q = 2J d , g 1(2) = (J z ^ ab =f Jl, &)/2 Jo- 
Fig. 3 presents a dispersion of the spin excitations in even and odd channels for both 
configurations of a spin alignment for model with parameters g\ = 0.05, g 2 = 0.15. 
Consider them in light of the above suggestion about simultaneous transformation F z — > 
AF Z of the spin configuration and superconducting transition in bilayer cuprates. If it is 
true for our spin model (3), then at T < T c the intense peak appears in the odd channel of 
X° dd at q ~ (n, n) with the frequency u> = A odd = 2fi v /gi, equal to a gap in the excitation 
spectrum in this channel for AF Z spin configuration. This peak is absent in case of F z 
configuration which is supposed to be ground state configuration at T > T c . Note that for 
F z configuration the low frequency excitations in x"(g,c<j — > 0) correspond to momenta 
q distributed along the circle |q — Q| = \f T 2g[ instead of the discrete incommensurate 
momenta Q v = (7r(l+r]), ir), observed for the low frequency excitations in YBCO |T3]].They 
imply an existence of the incommensurate spin structures. Fig. 3c schematically presents 
the branches of the spin excitations for the spiral states or the states with modulation 



of a local spin < S zn >= rfocosQ^n of the monolayer (1L) model [p0| . In |20| the 
resonance frequency have been identified as a frequency at the crossing point of branches 
at q = (n, 7r). 

Thus, it is confirmed that local magnitization do is retained in large doping range for 
wide diapason of parameters t', t z . The crossover of two spin configurations of the bilayer 
system proceeds at some doping. The lower normal state is that with the AF Z or F z spin 
alignment at small or large doping correspondingly. Maximum of the energy difference of 
these configurations corresponds to the optimal doping of the parent monolayer system 
and is connected with the maximum splitting of van-Hove singularities. The bilayer 
splitting is zero or small for the AF Z configuration, but large for the F z spin alignment. 
The observed in BSCCO large bilayer splitting in normal state, but small splitting in 
SC state H] may indicate according our models that the transformation F z — > AF Z of 
configuration proceeds simutaneously with supercunducting transition. It is shown that 
the magnetic excitation spectrum depends dramatically on type of a spin configuration. 
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The model spin system is presented for which a change of configuration F z — > AF Z is 
accompanied by appearance of resonance (or spin gap) in x° dd at g = (7T, 7r). 

A check of the above hypothesis requires further study, in particular calculations of 
magnetic excitations in case of the inhomogeneous (of the stripe type) spin structures of 
bilayer systems. 

Work is supported by Russian Fund of Fundamental Reasearch (Projects No. 00-03- 
32981 and No. 00-15-97334. Athours thanks V.Ya Krivnov for useful discussions and 
P.Fulde for possibility to work in Max Planck Institute for Physics of Complex Systems, 
Dresden. 
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Captions to Figures 

Fig. 1. Doping dependence of the local magnitization d (S) (top), of the T C (S) (middle) 
and value Wi(5) of the anamalous order parameter (bottom) for configurations AF Z F z 
(curves 1 and 2 correspondingly). Curve 3 is the difference of energies of these configura- 
tions. All values T c , t z , t', AH are in unit t. 

Fig. 2. Density of state for two types AF Z and F z of spin alignment for models with 
small t z ,t'. Thin solid and dashed lines refer to contributions in DOS from bonding 
and antibonding bands. The calculated DOS is smoothed with gaussian function with 
dispersion 8E = 0.02t, E is in unit t. 

Fig. 3. Changes in the magnetic excitations branches ui(q) during transformation 
F z — > AF Z of spin configuration of model (5) with anisotropic interlayer interaction with 
parameters g\ = 0.05, g-i = 0.15. A right figure is schematic presentation of the spin 
excitations in case of the spiral or modulated incommesurate spin structure with Q v = 
(7r(l — 77), 7r), 77 = 0.1. Squares mark the supposed resonance frequencies in various 
interpretations. 
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